AeroChem TN - 202 


DOE /JPL 95ho 62 - 79 /5 
b stribution Category UC-6i 


DEVaOPiv’:£Nr Of A MODEL 
AND COMPLETER CODE TO 
DESCRIBE SO* GRA*)E SiilCON 


PRODUCTION PROCESSES 


FIFTH QUARTERLY REPORT 
R. SRIVASTAVA AND R.K GOULD 
FEBRUARY 1979 



The JPL Low-Cost Solar Array Project is sponsored bv the 
U. S. Department of Energy and forms part of the Solar 
Photovoltaic Conversion Program to initiate a major effort 
toward the development of low-cost solar arrays. This work 
was performed for the w^et Propulsion Laboratory, California 
Institute of Technology by agreement between NASA and DoK. 

(NASA-C8-158601) DEVEtOFnENT OF A l^CDEl AND N79-23‘^14 

COHPOTEB CODE TC DESCRIBE SOLAR GRACE 
SILICON PRODOCTICK PROCESSES Cuacterly 

Report (AeroCheni Research Labs., Inc.) 36 p Onclas 

HC A03/HE AOl CSCL 10A G3/44 25172 


Research Laboratories. Inc. 



Princeton, New Jersey 



This report vas prepared ae an account ot work sponsored by the United States 
Government. Neither the United States nor the United States Department of 
Energy, nor any of their employees, nor any of their contractors, subcontractors, 
or their employees, makes any warranty, express or Implied, or assumes any legal 
liability or responsibility for the accuracy, completeness or usefulness of any 
Informatloii, apparatus, product or process disclosed, or represents that its 
use would not infringe privately owned rig'^ts. 



AeroChem TN - ZOZ 


DOE/JPL 954862 
Distr.bu'tion Cofegory UL-63 


KVELOPMENT OF A MODEL 
AND COMPUTER CODE TO 
DESCRIBE SOLAR GRADE SILICON 
PRODUCTION PROCESSES 


FIFTH quarterly REPORT 
R. SRIVASTAVA AND RK. GOULD 
FEBRUARY 1979 


JPL Confrocf Na 954862 
DRL Item No. 4, DRD No.QR 
LSA Silicon Moteriol Task 


Ap proved 











Hortwell F Coicote 
Director of Reseorch 



Research Laboratories, Inc. 


Princeton, New Jersey 



FOREWORD AND ACKNOWLEDGMENTS 


This report covers the period 1 October 1978 to 31 December 19" ~ Tu.rir’ 
this period, D.B. Olson and P.J, Howard contributed to the program. 

One of che authors (RS) wishes to express his gratitude to Professor D.E 
Rosner of Yale iniversity, for many helpful discussions pen aining ^ 
properties, and to R. Israel (Yale University) for providing Chopman-Er* Vog 
theory estimates of a . 



TN-202 


ABSTRACT 

This program aims at developing mathematical models, and computer codes 
based on these models, which will allow prediction of the product distribution 
in chemical reactors in which gaseous silicon compounds are converted to conderiSed- 
phase silicon. The reactors to be modeled are flow reactors in which silane or 
one of the halogenated silanes is thermally decomposed or reacted with an alkali 
metal, Hz or H atoms. Because the product of interest is particulate silicon, 
processes which must be modeled, in addition to mixing and reaction of gas-phase 
reactants, include the nucleation and growth of condensed Si via coagulation, 
condensation, and heterogeneous reaction. 

During this report period computer codes were develo: ed and used to cal- 
culate (i) coefficients for Si vapor and Si particles describing transport due 
CO concentration and temperature gradients (i.e.. Pick and Soret diffusion, 
respectively), and (ii) estimates of thermochemical properties of Si n-mers. 

The former are needed to allow the mass flux of Si to reactor walls to be cal- 
culated. Because of the extremely large temperature gradients that exist in 
some of the reactors to be used in producing Si (particularly the Westinghouse 
reactor), it was found that thermal (Soret) diffusion can be the dominant trans- 
port mechanism for certain sizes of Si particles. The therraochemical estimates 
are required to allow computation of the formation rate of Si droplets. With 
the completion of these calculations the information and coding of the particle 
routines in the modified LAPP code is at the point where debugging can be done 
and that is now in progress. 
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I, INTRODUCTION 


This program aims at developing mathematical models, and computer codes 
based on these models, which will allow prediction of the product distribution, 
as a function of time, in chemical reactors in which gaseous silicon compounds 
are converted to condensed-phase silicon. The reactors to be modeled are flow 
reactors in which silane or one of the halogenated silanes is thermally decom- 
posed or reacted with alkali metals. Because the product of major interest is 
condensed-phase silicon, the models considered must allow calculation of par- 
ticle formatior- (via chemical reaction followed by nucleation), and the growth 
rate of the particulate Si via coagulation, condensation, and heterogeneous 
reaction. Since the efficiency of the particulate Si collection process will 
depend on the particle size distribution of the Si, this must also be computed. 

During this quarter the following problems have received attention: (i) 

continued development of computational schemes for calculating heat and mass 
fluxes to walls of reactors in which SiXi./M (X = halide, M = alkali metal) 
flames are occurring, (ii) development of nucleation models by which Si(i) 
particles are formed from Si(g), (iii) continued development of detailed chemi- 
cal kinetic mechanisms describing the formation of free Si (Si(g)) by these 
flames, and (iv) development of coagulation and condensation models by which 
the Si(i^) particle growth is calculated. The bulk of this report concerns the 
development of models used to obtain the various transport coefficients needed 
to allow Si mass fluxes (as Si(g) and as Si{l) droplets) to reactor walls- 
Particular attention is paid to the role of thermal (Soret) diffusion which, 
because of the extremely large temperature gradients and the size of Si drop- 
lets that may be encountered in these reactors, needs investigation. Finally, 
the methods used to estimate the thermochemical properties of Si n-mers, 
needed to implement the Si nucleation model, are outlined. Short discussions 
of work on (iii) and (iv) above are also included. 

II, SILICON TRANSPORT LAWS AND IMPLICATIONS 


The mechanism of deposition of silicon particles onto reactor walls is 
dependent on the size of the particles. For particles larger than about a 
micron, deposition is controlled by turbulence-induced inertial impaction; 
sub-micron sized silicon droplets deposit by a convect ive-dif fusion process. 
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In order to predict the mass flux to the reactor wall of the latter class of 
droplets it is necessary to recognize two separate contributions to the time- 
averaged diffusional flux: (i) a concentration gradient-driven (Pick) flux and 

(ii) a temperature gradient-driven (Soret) flux* The relative importance of these 
contributions itself varies with particle size, the former being dominant for the 
smaller sub-micron sizes while the latter controls the behavior of particles with 
sizes closer to a micron, as shown qualitatively in Fig. 1. Thus, in general, 
for a species i (vapor or particle) that is dilute in a carrier gas the dif- 
fusional mass flux (per unit area, per unit time) at any point, jV, (correspond- 
ing to a particular pressure and temperature) is given by^*^: 


y: » - CD. |(TY.) + a.Y. (1 - Y.) J (7T) 

1 ill i 1 1 T 


turb 


( 1 ) 


where p is the local mass density of the carrier gas, the local mass fraction 
of the diffusion species, T the local absolute temperature (with the operator V 
denoting gradient of the particular scalar, temperature or concentration). The 
first two terms in Eq. (1) represent the diffusional mass fluxes due to concen- 
tration and temperature gradients, respectively. The third term, is the 

mass flux contribution due to turbulent or **eddy diffusion.*'^ Assuming, as is 
customary, that the turbulent diffusion also has a Fickean nature, one can write 


turb 


-DC (7Y.) 
P 1 


( 2 ) 


where e is the so-called **eddy dif fusivity Particles larger than about a 
P 

micron would inevitably be deposited on the walls due to inertial impaction 
caused by turbulent velocity fluctuations normal to the wall. Hence, for such 
large particles the overall mass flux to the wall would be dominated by 
However, the sub-micron size particles will b^'. driven close to the wall by tur- 
bulence eddies before the concentration and temperature gradient driven fluxes 
(i.e., Fick and Soret effects) become important in the final stages of deposi- 
tion, within a nearly laminar (viscous) sublayer next to the wall. 

However, before deposition rate (or mass flux) can be calculated, one 

needs to establish the Fick (or Brownian) diffusion coefficient D., the dimen- 

1 

sionless tiiermal (Soret) diffusion or *’thernophoret ic" factor and the "eddv” 
diffusivlty The former two transport property coefficients depend upon the 

molecular nature of the species present in a give, ’’effective” binary mixture 
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(e.g., silicon-argon, in the case considered later in this report), and on the 
prevailing pressure and temperature.* In addition, they are both strong func- 
tions of particle size, changing by several orders of magnitude (as shown later) 

+ 

in traversing the particle size spectrum from molecular diameters' to about a 
micron. 

Interestingly enough, while the values of the Pick dif fusivities (i.e., 
and e^) are always positive, ensuring that Brownian and turbulent transport 
of Si always occurs in the direction of decreasing concentration (towards the 
reactor walls), the value of may sometimes be negative. The latter circum- 
stance is usually characteristic of diffusing species that are lighter than 
the carrier gas (e.g., silicon vapor (28.09) diffusing through argon (39,94)). 
Thus the Soret flux has the potential to drive silicon vapor towards higher 
temperatures, away from the “cool** walls, and thereby reduce the overall col- 
lection efficiency of Si. ^ Fortunately, as shown later, the value of for 
Si particles (defined here as any molecular cluster exceeding roughly twice 
Che mass of a single vapor molecule) is invariably positive, since the parti- 
cles are heavier than the surrounding gas. That is, while condensed Si drop- 
lets w^ill always be driven towards the reactor walls, uncondensed Si 
could actually be driven away from the walls due to the Soret effect. In this 
report, therefore, in addition to considering the mechanisms of Brownian and 
turbulent transport of silicon, we investigate the conditions under v’hich Soret, 
or temperature gradient-driven diffusion, can be important. Although this 
latter effect has been ignored in the past in the belief that its contribution 
to the overall mass flux would be negligible, we conclude that the Soret effect 
can indeed play a crucial role in the Si (particle/vapor) separation process, 
under typical reactor conditionL.. Fortuitously, the Soret (or thermophoretic) 
flux for condensed Si droplets will always favor deposition onto the walls. 

The detrimental tendency of this flux to oppose deposition of any uncondensed 
Si vapor can be minimized, as sho\ later, by enriching the reactor with a 
light gas like H 2 , which tends to drive the system towards a condition of 
positive "effective thermal diffusion." 


A 

Under typical silicon re/'ctor operating conditions the concentration 
dependence of these transport coefficients will be negligible, cue to the 
diluteness of silicon in the carrier gas. 

Particle size is described in terms of diameter or radiftis since it will 
be assumed here that the particles (silicon liquid droplets) are perfectly 
spherical . 
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III. DESCRIPTION OF THE TRANSPORT COEFFICIENTS 


For the sake of comparison, it is necessary to define the ”eddy*' diffusivity, 

e , for mass transport. In general, by viewing turbulence as the random agitated 
P 

movement of small packets of fluid, it is possible to define c in terms of a 

P 

‘‘mixing length" (in much the same spirit as the ‘‘viscous mean-free-path" in the 
kinetic theory of gases represents an average distance for effective transfer 
of momentum by molecules). Thus, 


£ u'iil ici 


liHi 

•dy' 


where £i is a mixing length, u' the large scale turbulent velocity fluctua- 
tions present in the core flow of the reactor, and the magnitude of the 

dU 

time-averaged (steady) normal velocity gradient. Both i\ and |^j will vaiy 
with distance above a wall, y, with being zero at the wall since turbulence 
is precluded (i.e., the "no-slip" condition). While the many available models 
of turbulence provide different expressions for the variation of li with dis- 
tance normal to a wall, accurate estimates of and the proportionality constant 
in the above expression, for any given flow situation, usually result from ex- 
periments. Lin et al^ suggest the following expression, valid in the viscous 
sublayer flow, next to a wall: 


(ib)' 


< 5 


(3) 


where v is the kinematic viscosity and y is a non-dimensional normal distance 
from the wall defined as: 


y^ "1^ U(f/2)^^^ /v (4) 

with U the average velocity and f the Fanning friction factor. For particles 
small enough to follow the small-scale eddy motions close to a wall, the eddy 
diffusivxty does not depend upon particle size. 

Although silicon particles of widely varying sizes must be considered in 
the modeling of the deposition processes within any proposed reactor, reliable 
information about the "molecular" transport coefficients (i.e., and a^) is 
available only at the two extremes of the size spectrum. For particle si*» 
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much smaller than the mean-f ree-path of the carrier gas, the well-known 
kinetic theory results of Chapman and Enskog* (CE) apply. In this limit of 
**free-molecu3 e*' flew, surrounding gas molecules do not suffer appreciable 
changes iv ■’Vtribution functions upon collision with the particle. 

Firtherm since <pressions for are an outcome of CE theory using 
second-order terms of the perturbation series, while the expressions come 
out of first-order terms, one concludes that will be more sensitive to 
interactions between colliding molecules. Thus more sophisticated molecular 
potential models than the Lennard-Jones 12:6 (e.g., the exponential model) 
might ultimately be necessary to accurately determine in this limit, al- 
though the Lennard-Jones might suffice for engineering purposes. 

For particles much larger tnan the raean-f ree-path of the surrounding gas, 
the gas structure becomes "invisible” and the behavior can be described by the 
Navier-Stokes equations for a "continuum," subject to the usual "no-slip" 
boundary conditions at the particle surface (e.g,, the Stokes-Einstein (SE) 
expression for particle dif fusivity^) . 

The main difficulty in specifying transport properties is encountered 
for particles in the intermediate size range, on the order of a mean-free- 
path. This presents the notorious unsolved problem of the "transition regime." 
Although several attempts have been made to extend the continuum results into 
the transition regime, by imposing "slip" or jump-type boundary conditions at 
the particle surface, these are at best h»^uristic. Understanding of the real 
transition mechanisms involved is altogether insufficient. However, the transi- 
tion regime experimental data of Millikan and others have made it possible to 
estebiish various empirical correction factors (as discussed later) that cover 
a wide portion of the transition regime. Unfortunately, even such widely used 
interpolation formulae as the Stokes-Einstein-Millikan (SEM) expression are 
unsatisfactory since they do not exactly match the predictions of I'E theory 
when extended to the "f ree-raolecule" limit. 

Crucial to the accurate modeling of silicon deposition, therefore, is 
a proper description of the abovement ioned transport coefficients. Since 
present knowledge regarding the "transition regime" is insufficient, it would 
be justifiable to suggest universal interpolaticr. formulae for t:\ese coefficients. 
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covering the entire particle size (or Knudsen number*) range. Moreover, 
mathematically continuous functions that blend together the known results in 
different regions would be both physically and computationally more desirable. 

In what follows, a novel rational approach^ to creating such "universal” 
descriptions of the Pick diffusivity (i.e., D^) and the Sorer factor (i.e., 
a^) is outlined. Finally, some numerical results are presented to demonstrate 
the validity of the proposed formulae under conditions typical of a silicon 
reactor. ^ 

A. UNIVERSAL FORMULA FOR 

The expression for diffusivity adopted in t\is report blends the near^ 
continuum SEM relation (valid for small Kn^ values) with the f ree-molecule CE 
expression (valid for large values of Kn^, > 10) .according to the following 
universally valid relationship: 

^ ^CE * ^SEM 

^i ^SEM [1 + exp {- C(Kn^ - Kn©)}] (5) 

where C " nc are constants. The magnitude of C controls the abruptnes^^ of 

the tra on from to as the Kn.- is increased, Kno is the Knudsen 

SE>! CE 1 

number value at the "point of inflection" of Figure 2 ill xstrares the 

rationale behind choosing such an interpolation function to describe :he cif- 
fusivity of particles in a gas. 

‘ ^toues-Einstein-Millikan expression for di ‘:sivity, based 

on M- -i s*° famous oil-drop experiments which permitted measurement of 
the i\ *Tial drag force on p irticles of various sizes moving through a gas. 
Using a stochastic analysis of Brownian notion, Einstein^ ^ had earlier related 
the particle diffusivity to the "mobility," or relative velocity per unit force. 


The particle Knudsen number (Kn^^) will be defined here as the ratio of the 
mean-f ree-path (O in the surrounding gas to the particle radius (i.e., 

Kni = • Thus Kn^ 0 corr€^r..nds to the "continuum" limit while Kn^ ^ 

refers to the "f ree-molecule" limit. It is some imes asrumed that the 
range 0 < Kn^ < 0.25 represents the region of "slip-flow" while 0.25 •' Kn^ ^ 
10 is the region of "transition-flow." 


6 
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through the well-known relation: 


= kTB (6)* 

the mobility^ B, being given by the Stokes formula for low Reynolds number 
continuum flow around a sphere: 

® ■ 6 ^ <” 

where ^ is the dynamic viscosity of the gas. Recognizing that as the particle 
size approaches the gas mean-free-path, the drag for a given velocity becomes 
less (due to ”slip**) than predicted by Stokes law (i.e., the mobility increases), 
Millikan suggested the following expression for mobilit\ , which is a well- 
accepted correction of the Stokes result: 


5 


1 


(7a) 


where 


£ H [1 + C Kn. + C,Kn. exc(- O.S8/Kn.)]“* (7b) 

m 1 G 1 ' X 

C = (2 - a )/a (7c) 

m mm 

(^) 8 "m' j 

with a denoting the so-called momentum accommodation coefficient, a measure 
m 

of the efficiency of momencuin transfer between the particle and the surrounding 
gas molecules. Combining Eqs. (^) through (7d) yields the needed expression 
for diffusivity, v^alid at least for small departures from Rn. = 0: 


"k 

k is the Boltzmann constant. 

Values of depend on t..^ nat re of the particle and the surrounding gas. 
Experimentally determined values for various substances^ usuallv lit* in the 
range 0.9 < < 1.0. In this report = 1 (perfect accommodation ) is 

assumed for che results '^^esented later. 
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D 


SEM 


kT 


[1 + C Kn. + C.Kn. exp(- 0.88/KnJ] 
m i a 1 1 


( 8 ) 


Ac the other extreme, for particles much smaller chan the gas mean-free- 
path (but not necessarily as small as molecules) the diffusivicy D-_ should be 
given by the Chapman-Enskog relation**: 


CE 




(9) 


where p is the pressure (atm) , M and the. molecular weights of the carrier 

gas, c, and diffusing molecular species, i, respectively, c . the effective 

ci 

collision diameter (angstroms) and H the reduced collision integral. Equation 
(9) is strictly valid only for molec’..:les with spherically symmetric potential 
force fields, tfich the supplementary relations: 


ci 


2 ^‘'c ®i^ 


n 



(See Ref. 13 for a tabulation 
of this function.) 


( 10 ) 


'ci 






where and can be taken to represent the size and energy«well depth in, 
say, Che Lennard-Jones molecular interaction model. 

To generalize Eq. (9) for the case of small particles, we need to redefine 
only the molecular weight and the mol-'-lar diameter Treating the par- 

ticle as a "heavy molecule," g times the mass of the corresponding molecule, 
its "molecular weight" can be expressed as: 



(ID* 


A 

N. is the Avogadro number, 

A 
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where Mi is the oolecular weight of the substance constituting the particle, 

g is the ratio of the particle volume to the molecular volume (i.e., ■ ) 

"fP i 

of the substance, and is the particle density. The "molecular" 
diameter of a particle is simply 

°i * ^^^i^ ** 


If the particle radius is expressed in centimeters. Combining Eqs. (9) 
through (12) yields the following expression for the diffusivity of very 
small particles: 



(13) 


where 


0 


2 

ci 


i (Cg + 2 X 10* Ri)* 


(13a) 


r. 



(13b) 


g 



(13c) 


Incorporating the expressions for and given by Eqs. (S) and (13), 

into Eq. (5) yields the required, universally valid, equation for the diffu- 
sivity The value of Kno was taken to be 0.44 so that the point of inflec- 

tion in Eq. (5) coincided with the inflection point of the SEM formula (i.e., 

Eq. (8)). The value of C was chosen to be 1.5 for the results presented belov. 


b. UNIVERSAL FORMULA FOR 


Using a rationale similar to the one adopted for D^, but necessarily 
different because entirely different (and more complicated) expressions 
usually describe we propose the following universally valid formula 
for the Soret factor: 


9 
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( ( 1 \ 
Vl + AKn, /“• ■*" \l + AKn, / 


Oo 


(14) 


where a» is the '’thermal diffusion factor" predicted by CE theory for Kn^ -► ® 
and Gio is the value of the Soret factor as Kn^ 0. The constant A is deter- 
mined by forcing the prediction of Eq* (14), in the transition regime, to 
agree with established experimental data for some value of Kn^ of order unity. 

It is known that a particle suspended in a gas with a uniform tempera- 
ture gradxent will drift through the gas under the action cf a thermal force. 
When this "thermophoretic" force just balances the opposing drag force on the 
particle, a constant thermal settling velocity is achieved. Using Millikan's 
oil-drop apparatus, modified to provide for a temperature gradient, Phillips^** 
was able to accurately describe the thermal settling velocity up to » 0(1), 
for particles of widely varying thermal conductivities (including the coupling 
that usually exists between the velocity and temperature fields surrounding a 
moving particle). We have incorporated^ the Phillips formulae into a computer 
program that predicts the thermal settling velocity (V^) at some value of the 
Knudsen number (say, Kn^ « 1) and used this to find the corresponding Soret 
factor (ap). Knowing ajo and a©» as shown below, A can then be determined as: 


A 


1 — Jt o 


(15) 


for a given particle-carrier gas combination and given flow conditions. 

Since the CE expression for is extremely complicated, we utilize the 
analysis of Waldmann and Schmitt^’ which yields the following expression for 
valid for small particles whose drift in a temperature gradient can be 
associated with the stronger bombardment by higher energy molecules from the 
"hot side": 


1 

T 


(VT) 


Y H a.D. = f ) 

1 ‘ 


D. 
« 1 


where is the diffusivity already discussed, v the kinematic viscosity, and 

a, the fraction of molecules that undergo "diffused reflection" from the par- 
a 

tide surface. It turns out that the above result for a. does not exactlv 

1 
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match the CE result for molecular sized particles- In order to achieve this 

match a new constant correction factor o is introduced as follows: 

corr 


a 



1 

a 

corr 


(16) 


Physically, tha effect of account fcr the variation in with 

increasing particle size.* The rat'^o (v/D^'^ is che particle Schmidt number 
(Sc^). It is a measure of the relative importance of momentum to mass trans* 
fer. Since the diffusivity of large particles is small, unlike gas molecules, 
Sc^ is typically large (as seen later. Table I). 

Here olq is taken to be related to the thermc horetic velocity deduced 
by Epstein,^’ for particles much larger than the .:as mean-free-path (using 
a slip-flow procedure). Using Epstein*s result oAe can write: 

3 . [ 1 

- 4 i I 1 + 0.5(X./A ) 

* 1C 

where and are the thermal conductivities of the particle and gas, 
respectively. 

Equations (15), (16), and (17) when incorj orated into Eq. (14) provide 
the required universal formula for the Soret factor 



C, TRANSPORT PROPERTIES 


1. General Trends 

The analysis just presented for the evaluation of the transport 
properties of particles dilute in a carrier gas, was first applied to the 
case of silicon particles in argon, under typical reactor operating condi- 
tions. Figures 3 and 4 show the numerically computed values oi and 
over a wide range of Knudsen number (Kn^). Table I is p v-.ded to gain an 
idea of the corresponding particle size in terms of i other possible size 
variables: g, R^, or Sc^. Note that while g and are absolute measures 


* 


Vv-hen 


is negative, 1 may be assumed. 
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of particle size, Kn^ and Sc^ depend strongly upon the prevailing pressure 
and temperature. Table II gives the values of input constants and the com- 
puted values of some useful quantities. The following features regarding 
Figs. 3 and 4 are also noteworthy: 

(i) Both D and a. change rapidly by orders of magnitude, in going 
1 i o 

from very small (g ■ 1, ■ 1.58 A, Kn^ ■ 4.4 x lO , SCj^ = 0.73) to large 

particles (g«10^*, = 157 um, Kn^ = 4.4 x 10"’, Sc^^ - 4.6 x lO*). In this 

connection it should be mentioned that particles of about 1 pm diam (see Fig. 
1 and Table I), which lie within the "transition regime," correspond to the 
size ranges: 

10*® < g <10** 

X 10"’ cm < < 7.3 X 10”’ cm 

2.04 > Kn^ > 0.95 
1.1 X 10* < Sc^ < 5.35 X 10‘ 

(11) As expected physically, the dlffuslvlty of large particles Is 
small. Increasing significantly as particles approach molecular dimensions. 

On the other hand, the magnitude of the Soret factor shows an opposite 
trend In this case. The physical Implication^ of this Is better understood 
from £q. (1), which gives the net mass flux* Since is always positive a 

negative (or temperature gradient) tends to reduce the flux* In this 
case, the usually large when Is small, and vice versa, leads to the 
conclusion that the Soret flux will dominate the Brovmlan flux only for 
larger particles. 

2. Pressure and Temperature Dependenc es 

Figure 5 reveals the pressure and temperature dependences of the 
transport coefficients. Due to the complicated functional dependence of 
and on Kn^ and simple representations may not be possible for all 

particle sizes. It is clear, however, that over the range of temperatures 
( *^ 1700 to 3500 K) and pressures (» 0.5 to 1.0 atm) of Interest In silicon 
rc ors, retains approximately the same characteristics as described by 
I. £ expression (i.e., T^^*/r), at least for particle diameters up to about 
1 pm. On the other hand, Is nearly Insensitive to temperature and pres- 
sure. Figure 5 Is restricted to particle diameters up to about a micron 
since it is generally true that the Soret effect (thermophoresis) may not 
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be as significant as inertial forces for very large particles* 

3. Need for Effective Transport Properties 

Figure 6 compares the case of a negative (e*g*. Si vapor in Ar) 
versus one with a positive (e.g.. Si vapor in H 2 ).* In order to increase 
the separation (or collection) efficiency of silicon, it is clear that 
should be positive so that the Soret flux is directed towards the "cold** 
reactor walls* Since a multicomponent mixture of gases is actually present 
an ''effective** Soret factor must be considered such that 

“l.eff ■ [ Vl ]/“ ■ "l’ 

where is the mole fraction of species i and the Soret factor for species 
i diffusing in a carrier gas j. It is to be noted that the "effective" Soret 
factor depends both on the relative concentrations of the various gaseous 
species present as well as on the magnitude and sign of their individual 
Soret factors* Typically, as in the Westinghouse process, for instance, 
one can expect a mixture of silicon in mainly sodium chloride vapor, argon, 
and hydrogen gases* The various mole fractions would be roughly 


*Si 

= 0.08 


= 0.5 

*Ar 

= 0.12 

^aCl 

« 0.3 


If the effective carrier gas is a mixture of Ha* Ar, and NaCl its average 
molecular weight would be 


M 


avg 



23.34 


(19) 


Since the molecular weight of Si vapor molecules is larger than this (M^^ 28) 

one would expect a positive value of to govern the Soret diffusion of 


* 

Note, 


the value of for Si particle s in both Ar and H 2 is positive. 
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silicon vapor In a mixture of the abovementloned gases. Indeed, using the 
above mole-fraction estimates together with the following approximate values 
of a (deduced from C£ theory) applicable to Si vapor: 


“SI.H, 

“si.NaCl 

®Si,Ar 


0.7 


-0.3 


- 0.1 


one obtains 0.27. Thus the net Soret flux of vapor should normally 

be directed towards the wall. Using this it was estimated that the Soret 
mass flux for Si vapor to the reactor wall may not exceed 15-20% of the cor- 
responding Pick flux. However, the sign of a is sensitive to the pre- 

et t 

vailing product concentrations. It is conceivable, therefore, that in certain 
portions, near the reactor walls, the concentration of H 2 may be locally small 
enough for to become negative and lead to a detrimental flux of Si 

vapor away from the walls.* To eliminate the occurrence of such a possibility 
we would recommend maintaining a Ha rich atmosphere near the reactor walls 
at all times. 

Finally, it is worth pointing out that the real problem, as mentioned 
earlier, involves (besides vapor) some silicon droplets that may be treated 
as "heavy molecules" (i.e., the submicron size class). In the context of the 
above discussion, pertaining to vapor, it is clear that the effective "molec- 
ular weight" of this "heavy-molecule-gas" will always be considerably greater 
than that of Si vapor (see £q. (11)). Hence, one can, in fact, expect a huge 
(but favorable) Soret contribu*. -un to the overall mass flux for particles in 
the "transition regime." 

4. Relative Importance of Various Transport Mechanisms 

Finally, one needs to gain some idea of the relative importance of 
the three mechanisms of diffusion (i.e.. Brownian, Soret, and turbulent) 
studied in this report* Brownian and turbulent diffusion are known to be 
important for very small and very large particles, respectively. However, 


Our more exact boundary layer calculations in the future will also be 
aimed at determining the extent of this possibility. 
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the importance of Soret diffusion for particles of about a micron diameter is 
easily seen from Fig. 5. Typically, is large ( ■■ 10^) while is small 
( •» 10"*). Recalling that the Soret factor, is related to a thermal dif- 
fusion coefficient D_ . (which is really the appropriate counterpart of D ), 

* » 1 X 

the following useful ''thermal diffusion ratios" may be defined: 

^ ^ “ “I’^i « - ''i> <“) 


4 




p(Di + Ep) 



"i°i^i 



( 21 ) 


The former compares the extent of Soret and Brovmian diffusion, whereas the 
latter expresses the importance of Soret vs overall Fickean (i*e«, Brotmian 
turbulent) diffusion. Since, as seen above, a.D. » 1, Y. a 0.1, and e 10“* 

111 p 

(from £q. (3)), close to the wall,* one finds: 


« 10 " 

k^ 10 

That is, Soret transport of silicon particles will predominate all the other 
forms of diffusion for particles about 1 pm diara, in a viscous sublayer region 
close to the reactor wall. 

5. Particle Size Range Over Which Soret Transfer Dominates 

It is of interest to determine here the particle size "window" (see 
Fig. 1) for which Soret diffusion dominates. An upper limit is approximately 
set by the physical requirement that particles exceeding about 1 pm diam will 
be controlled solely by their inertia.* However, the condition k,p ® 1 (see 
Eq. (20)) determines the lower size limit at which Soret diffusion becomes as 
important as Brownian diffusion. It is interesting to note that this lower 
limit is sensitive to the prevailing Si concentration (or mass fraction). 


Note that since Cp w y^, the value of Cp drops rapidly as the wall is 
approached. For 0 < y“*" < 5, the variation in Cp corresponds to 0 < Cp < 
0.15 (see Table II). 
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Using obtains which (from Fig. 5) corresponds to a 

particle diameter of 0.01 um» under typical reactor conditions. With » 10 
one can expect Soret diffusion to outweigh Brownian diffusion, a condition 
which here would correspond to » 100 and a particle diam of 0.02 um. Thus 
one can safely say that Soret diffusion (i.e., thermophoresis) will be the 
dominant transport mechanism, near the reactor walls, for silicon droplets in 
the size range: 0.01 ym to 1.0 pm diam. Note that this is a much wider size 

range than the one usually quoted in the literature on a somewhat ad hoc basis. 

6. Conclusions 

The analysis just presented emphasizes the following aspects relevant 
to modeling the silicon separation processes within a reactor: 

(a) It is essential that the model describe the transport to the 
reactor walls of a whole size distribution of silicon particles, from vapor 
molecules to droplets of several microns in diameter. As a first step to 
determining the average mass flux due to a statistical particle size distribu- 
tion, the universally valid transport coefficient formulae, developed in this 
report, will be needed. 

(b) Very small or large silicon droplets will be deposited on the 
reactor walls by Brownian and turbulent diffusion, respectively. However, 
the often Ignored transport mechanism of Soret diffusion (or thermophoresis), 
due to large temperature gradients, will control the deposition of silicon 
droplets in the diameter range 0*01 pm to 1.0 pm, within a viscous sublayer 
region close to the walls. 

(c) Unlike the Brownian and turbulent diffusivities, the Soret 
factor can possibly take on negative values, rendering the Soret flux of 
silicon vapor to be directed away from the reactor walls. For silicon drop - 
lets , however, since this factor is inevitably positive, Soret diffusion 
invariably favors deposition onto the walls. It is to be noted that while 
the Soret flux of silicon vapor in a typical reactor may not exceed 15-20% 
of the prevailing Fick flux, close to the reactor walls, the contribution of 
the Soret effect to the overall mass flux of silicon droplets far outweighs 
all other diffusion mechanisms (as pointed out in (b) above). To eliminate 
the possible occurrence of a detrimental Soret flux of uncondensed silicon 
vapor in local regions close to the reactor walls, we recommend that a 
hydrogen-rich atmosphere be maintained within the reactor at all times. 
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IV. DESCRIPTION OF THERMODYNAMIC PROPERTIES 


A. SILICON N-MER THERMOCHEMISTRY 

In this report period we have also concentrated efforts on obtaining a 
set of thermodynamic data which will allow us to model tae formation of Si(£) 
from Si(g) in a reasonable fashion. The chemical kinetic model we are using^’ 
requires values of the heat capacity, c , the enthalpy (H^ - ** 298 ^’ 
energy function -(G® - H|jg)/T, and the heat of formation, AHf®, for each 
Si n-mer. Such data do not exist and must therefore be estimated. Two 
approaches have been chosen. The first uses a very naive chemical model to 
estimate the unknown quantities; it assumes that the heat capacity of each 
n-mer is simply the heat capacity of Si(i) times the number of Si atoms in 
the n-mer. The free energy function is obtained by assuming that the free 
energy of the n-mer is the sum of the molar free energy of the liquid (times 
n) and the molar surface energy, G^, where, from the Kelvin eq.,*“ 


s,n 


3 “ * 'n °SlV 


( 22 ) 


In this equation r^ is the radius of the n-mer and is assumed to be that of 
a liquid drop of Si containing n Si atoms with a density given by’® 

Og^ « 3.02 - 0.35 X 10 -^ (T - 273.15), gm ml’’ (23) 

The surface tension of Si liquid, Og^, is given by*® 

o = 720[1.67 - 0.67 (T/1685)’-*] ergs cm'* (24) 

wX 


No attempt is made to enforce self-consistency in the model, i^e., the entropy 
is not the integral of the heat capacity divided by T» The thermodynamic data 
in Table III are examples estimated using this model* Generally the model 
seems to suffer from underestimating the entropy (and thus overestimating free 
energy functions), making it less than zero (and thus not physically reason- 
able) at lover temperatures* However, this model does have the virtue of 
converging towards the thermodynamic properties of the liquid, as it should. 
Thus for large n-mers and high temperatures the model will be realistic* 
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The second model used to estimate n*mer thermodynamic data is that of 
Bauer and Frurip^^ also described in Refs* 17 and 23. This model can be ex*- 
pected to yield more realistic values of thermodynamic properties of n-mers 
where n ^ 100 and at low temperatures since the use of the bulk liquid sur- 
face energy is avoided* However » it does not guarantee convergence to the 
liquid thermodynamic properties at large n and, as can be seen in Table IV, 
the value of the free energy function actually moves away from the bulk liq- 
uid values as n increases* The strength of the model is obviously at low 
values of n where the classical model gives physically unrealistic results. 
Thus, for now, for n < 100 we will use the Bauer model estimates. For n > 

100 we switch to the classical model. These thermodynamic data, together 
with estimates of reaction rates obtained using the Bauer and Frurip kinetics 
model^^ are currently being used to compute Si n-mer growth in SiCl^/Na flames 
Since it is evident that universal formulae of the type discussed earlier 
in connection with the transport coefficients, will ultimately be needed to 
describe the thermodynamic properties as well, such efforts arc currently 
being pursued. 

B. OTHER WORK 

The last quarterly report^ ^ described the development of detailed mech- 
anisms to conserve computer run time. Further work during this quarter on 
simplifying the original 25 step reaction schemes^* for these flames has 
reduced the number of reaction steps necessary to describe these flames to 
13. Reactions involving minting of SiX^ species and H-atom stripping reac- 
tions (H + Six -► Six , + HX) do not appear to plav an important role in 
m m— 1 

the overall reaction scheme which yields free Si. 

The particle growth model is now inserted in the LAPP code*^ and 
debugging of this model has begun. 


V. PLANS 


During the next quarter, the SiCl^/Na reaction scheme, the Si n-mer 
thermochemical estimates, and the Si paiticle transport coefficients ob- 
tained in this quarter will be used to calculate heat and Si mass fluxes 
to walls in reactors like those being designed by West inghouse* ^ and 
AeroCl er.. ^ To do this we will be using the LAPP code,*’ now modified 
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to treat particle nucleation and coagulation in an enclosed reactor. The 
code is in the debugging stage. To obtain more detailed wall heat and mass 
flux estimates than are possible with the modified LAPP code we .^ill be rely- 
ing on the boundary layei code*’*’ which specifically examines the 
crucial region close to the wall. 


VI. NEW TECHNOLOGY 


No reportable items of new technology have been Ic'entified. 
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TABLE I 

OWPARISON OF VARIOUS SIZE RELATED VARIABLES FOR Sl-Ar 
(p « 1 atm, T » 2200 K) 


8 

Rj^ (cm) 

Kn 

i 


1.00 

1.58 (-08) 

4.39 (03) 

7.2 (-01) 

1.00 (01) 

3.39 (-08) 

2.04 (03) 

2.4 

1.00 (02) 

7.31 (-08) 

9.50 (02) 

7.8 

*.00 (03) 

1.57 (-07) 

4.41 (02) 

2 3 (01) 

1 00 (04) 

3.39 (-07) 

2.04 (02) 

1.1 (02 

1.00 (05) 

7.31 (-07) 

9.50 (01) 

5.2 (02) 

1.00 (06) 

1.57 (-06) 

4.41 (01) 

2.3 (03) 

1.00 :07) 

3.39 (-06) 

2.04 (01) 

1.0 (04) 

1.00 (08) 

7.31 (-06) 

9.50 

4.9 (04) 

1.00 (09) 

1.57 (-05) 

4.41 

2.3 (05) 

1.00 (10) 

3.39 (-05) 

2.04 

1.1 (06) 

1.00 (11) 

7.31 (-03) 

9.50 (-01) 

5.3 (06) 

1.00 (12) 

1.57 (-04) 

4.41 (-01) 

2.1 (07) 

1.00 (13) 

3.39 (-04) 

2.04 (-01) 

6.8 (07) 

1.00 (14) 

7.31 (-04) 

9.50 (-02) 

1.7 (08) 

1.00 (15) 

1.57 (-03) 

4.41 (-02) 

4.2 (08) 

1.00 (16) 

3.39 (-03) 

2.04 (-02) 

9.5 (08) 

1.00 (17) 

7.31 (-03) 

9.50 (-03) 

2.1 (09) 

1.00 (18) 

1.57 (-02) 

4.41 (-03) 

4.5 (09) 
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TABLE II 

VALUES OF SOME INPUT AND COMPUTED CONSTANTS 


INPUT CONSTANTS (Sl-Ar) 


P 

T 

“d 

Hx 

Me 

a„ 


“1 

(^l/k) 

(ec/k) 

C 

Kno 


) 


gm/gm-inol 


“ 1 atm 
» 2200 K 
= 1.0 
= 1.0 
= 28.09 
= 39.94 

■ -0.10279 (CE theory estimate) 
= 6.023 ^ 10** molecules /gm-mol 
= 1.3805 X 10““ ergs K“* 

= 2.910 A 
= 3036.0 K 
= 3.542 A 
= 93.3 K 
= 350 
= 1.5 
= 0.44 




Particle density: P 


mp 




i ^mp 
= 3.025 gm cm 
= 1.175 X 10" 
= 1685 K 


(See Eq. (5)) 

[1 - Z(T - Tnp)] 


Gas Viscosity: 


Ji 

Ur “ V^r/ 


Pr “ 8.021 X lO"*’ poise 

Tj. * 2000 K 

w *0.35 
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TABU II (Continued) 


COMPUTED CORRECTION FACTORS 


Cq 

Cd 

- 1.0 1 
» 0.61557 j 

> (see Eq. (8)) 

A 

> 1.9182 

(see Eq. (14)) 


COMPUTED VALUES OF Smffi VARTART.PS 


V 

= 3.7498 


p 

- 2.2116 

* 10“*’ gm cm”* 

£ 

= 6.9572 

X lO"* cm 

>c 

» 1.5472 

* lO”* cals/(cm-s-l 
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TABLE III 

THERMODYNAMIC PROPERTIES OF Si N-MERS 
USING A CLASSICAL LIQUID DROP MODEL^ 



c /n 

-(G® - H598)/nT 

(H^ - ^298^^^ 
kcal mol*^ 


cal mol“‘K“* 

cal mol”*K“* 

n 

(Si atOBj)”‘ 

(Si atom)”* 

(Si atom)*^ 


1000 K 


6 

6.5 


- 2.74 

4.56 

12 

6.5 


0.70 

4.56 

24 

6.5 


3.43 

4.56 

48 

6.5 


5.60 

4.56 

b 

6.5 


7.32 

4.56 

OB (liquid) 

6.5 

2000 K 

13.93 

4.56 

6 

6.5 


11.65 

11.06 

12 

6.5 


12.85 

11.06 

24 

6.5 


13.80 

11.06 

48 

6.5 


14.56 

11.06 

b 

6.5 


15.16 

11.06 

* (liquid) 

6.5 

3000 K 

17.47 

11.06 

6 

6.5 


17.70 

17.56 

12 

6.5 


18.13 

17.56 

24 

6.5 


18.47 

17.56 

48 

6.5 


18.74 

17.56 

b 

6.5 


18.96 

17.56 

•(liquid) 

6.5 


19.78 

17.56 


a 

AHf29s/^ for dll n-mers is 11.585 kcal mol (Si atom) 


b 


From Ref. 


21 . 
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TABLE IV 

THERMODYNAMIC PROPERTIES OF Si S-MERS 
USING THE MODEL OF REFERENCE 22 


n 

Co/n 

cal mi)l"*K“* 
(^^i fctom)”* 

-(G* - H|98)/nT 
cal mol-*ir» 
(Si atom)"* 

(Hj - H298)/n 

fecal mol"* 
(Si atom)"* 



1000 K 


6 

4.83 

16.97 

3.93 

12 

4.83 

16.90 

3.93 

24 

4.83 

17.64 

3.93 

48 

4.83 

18.72 

3.93 

96 

4.83 

19.88 

3.93 

• (liquid)® 

6.5 

13.93 

4.56 



2000 K 


6 

5.44 

19.89 

8.87 

12 

5.4i 

19.81 

8.87 

24 

5.44 

20.56 

8.87 

48 

5.44 

21.64 

3.87 


5.44 

22.80 

8.87 

X (liquid) 

6.5 

17.47 

11.06 



3000 K 


6 

5.28 

21.75 

14.31 

12 

5.28 

21.74 

14.31 

24 

5.28 

22.43 

14.31 

48 

5.28 

23.51 

14.31 

96 

5.28 

24.62 

14.31 

•(liquid)® 

6.5 

19.78 

17.56 


298.15 K 


AH|/n 

n Real Tnol~* (Si atoin)~* 


6 

73.00 

12 

63.23 

24 

55.01 

48 

48.09 

96 

42.29 

•(liquid) 

11.59 


^ From Ref. 21. 
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NUMBER OF MOLECULAR VOLUME UNITS, g 



FIGURE 5 VARIATION OF AND WITH PRESSURE 
AND TEMPERATURE 
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PARTICLE RADIUS (Ri), cm 


FIGURE 6 COMPARISON OF THE SORET FACTOR FOR 
TWO DIFFERENT TYPES OF CARRIER GASES 
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